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INTRODUCTION 


Atomic  oxygen  and  translational  temperature  play  vital 
roles  in  the  chemistry  and  heat  balance  of  the  earth's 
atmosphere  above  80  km  altitude.  However,  existing  methods  for 
measuring  thermospheric  profiles  of  atomic  oxygen  density  and 
translational  temperature  have  very  severe  limitations.  The 
paper  by  Zachor  and  Sharma  (1989),  an  expanded  version  of  the 
present  report,  discusses  the  chemical  processes  and  cooling 
mechanisms  that  involve  upper  atmospheric  atomic  oxygen.  It 
also  reviews  the  various  techniques  that  have  been  used  to 
measure  the  temperature  and  O-atom  density,  and  points  out 
their  deficiencies. 

The  present  study  was  undertaken  to  determine  the 

feasibility  of  recovering  vertical  profiles  of  translational 

temperature  and  oxygen  atom  density  from  measurements  of  the 

limb  radiance  profile  near  147  urn  and/or  63  uni  wavelength, 

.  3  3  3  3 

corresponding  to  the  01  transitions  (  Pq — *■  and  P^ — *•  P2) 

of  the  ground  electronic  state  of  atomic  oxygen.  The  assumption 
that  the  ^P  fine  structure  levels  are  in  Local  Thermodynamic 
Equilibrium  (LTE)  with  the  translational  temperature  is  crucial 
to  the  proposed  technique,  but  seems  a  reasonable  one  based  on 
the  very  long  radiative  lifetimes  of  the  level  transitions  (see 
Fig.  1).  The  principal  advantage  of  the  proposed  technique 
over  other  remote  methods  (see  Zachor  and  Sharma,  1989)  is  that 
no  models  are  needed  for  relating  the  observed  emission  to  the 
oxygen  ground-state  density. 
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The  sought-after  vertical  profiles  can,  in  principle,  be 
recovered  from  a)  a  pair  of  limb  radiance  profiles  representing 
the  total  (spectrally  integrated)  apparent  intensities  of  the 
two  01  lines  (Sharma,  et  al.,  1988)  or  b)  a  spectrally  resolved 
limb  radiance  profile  for  just  one  of  the  lines.  The  use  of 
spectrally  resolved  data  implies  an  instrument  with  resolving 
power  greater  than  2  x  10"*.  This  value  corresponds  to  resolved 
spectral  intervals  at  63  pm  and  147  pm  that  are  comparable  to 
the  respective  line  widths.  We  estimate  that  a  state-of-the- 
art  Fabry-Perot  interferometer  based  on  metal  mesh  etalon 
reflectors  could  achieve  resolving  power  as  high  as  5  x  10~*  or 
possibly  10^  at  these  wavelengths.  An  IR  heterodyne  system,  if 
one  could  be  designed  to  operate  near  147  pm  wavelength,  would 
be  more  complex,  but  would  easily  provide  the  required  high 
resolving  power. 

The  proposal  to  remotely  sense  both  temperature  and  0-atom 
density  using  a  spaceborne  Fabry-Perot  sensor  system  with 
63  pm  and  147  pm  channels  was  originally  put  forth  by  Sharma, 
et  al  .  ,  (1987,  1988).  The  1987  paper  outlines  a  preliminary 
sensor  system  concept  and  describes  its  critical  components. 
Sharma,  et  al . ,  (1988)  describe  the  'non-resolved  approach"  and 

a  method  for  inverting  the  dual  radiance  profiles,  and  show 
temperature  and  0-atom  density  profiles  retrieved  from  noise- 
free  synthetic  data.  We  have  since  developed  an  improved 
retrieval  algorithm,  applied  it  to  noise-contaminated  synthetic 
data,  and  on  the  basis  of  the  results  have  evaluated  the 
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feasibility  of  this  approach.  The  more  recent  work  will  be 
reported  in  the  forthcoming  Part  2. 

The  present  report  concerns  the  "spectrally  resolved 
approach".  We  describe  a  straightforward  retrieval  methodology 
and  the  characteristics  of  solutions  obtained  from  synthesized 
limb  spectral  data  contaminated  by  noise.  Included  in  the 
results  is  a  nomograph  that  defines  the  combinations  of 
instrument  noise-equivalent  spectral  radiance  and  resolving 
power  for  which  the  O-atom  density  and  translational 
temperature  can  be  successfully  recovered.  These  basic  data 
requirements  can  be  translated  by  the  experimentalist  to  design 
tradeoffs  and  system  parameters  in  order  to  gauge  the 
feasibility  of  a  given  sensor  concept.  In  this  fashion  we 
evaluate  the  feasibility  of  the  proposed  (Sharma,  et  al  .  ,  1987) 

Fabry-Perot  experiment. 

We  begin  with  a  summary  of  the  governing  radiative 
transfer  equations  and  a  general  description  of  the  retrieval 
method.  Section  3  gives  some  details  of  the  retrieval  algorithm. 
The  remaining  sections  present  results  and  conclusions. 
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2.  RADIATIVE  TRANSFER  EQUATIONS  AND  RETRIEVAL  METHODOLOGY 


The  transfer  of  photon  spectral  radiance  Nu  arising  from  a 
single  transition  at  wavenumber  vQ  through  a  slab  of  thickness 
Az  along  the  line-of-sight  can  be  evaluated  by  the  equation, 


Nv(z+Az)  =  Nu(z)  exp  (-Atu)  +  [2c,U02y/ (1-y)  ]  [1-exp  (-Axv)  ] ,  (1) 


where  c  is  the  speed  of  light.  Equation  (1),  for 
appropriately  defined  y,  is  not  restricted  to  the  LTE  case; 
see  Zachor  and  Sharma  (1985) .  In  the  present  case  LTE  is 
assumed,  and  y  reduces  to 

y  =  exp <-c2u0/T)  ,  (2) 

where  c2  =  1.4388  K/cm  1  is  the  second  radiation  constant 
and  T  is  the  slab  temperature.  The  first  factor  of  the 
second  term  in  the  r.h.s.  of  equation  (1)  is  then  equal  to 
Planck's  function  evaluated  at  u0,T,  while  the  second  factor 
represents  the  fraction  of  radiation  escaping  the  slab.  The 
monochromatic  optical  thickness  of  the  slab  is 

Axv  =  n  Az  f(V-\>0)  S(\)0,T),  (3) 

where  n  is  the  species  (O-atom)  density,  f (u)  is  the 
normalized  Voigt  shape  of  the  line  due  to  Doppler  and 
collisional  broadening,  and  S(v0,T)  is  the  integrated  intensity 
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of  t  3  line,  given  by 


S  (U0,  T)  =  (A/8JCC\>02)  [  (1-Y)  /y]gu  exp{-C2Eu/T)  /Q(T)  .  (4) 

Here,  A  and  Eu  are  the  Einstein  A-coef ficient  (s_1)  and  upper 
state  energy  (in  cm-1)  ,  respectively;  they  have  one  of  the 
values  given  in  Fig.  1,  depending  on  whether  u0  =  68  cm  1  (X0  = 
147  nm)  or  =  158.5  cm'1  (XQ=  63  ^m) .  The  statistical  weight 

gu  of  the  upper  state  is  unity  for  the  147  Jim  line  and  three 

for  the  63  |im  line.  The  partition  sum  Q  (T)  is 

Q  (T)  »  exp(-c2  x  226. 5/T)  +  3  exp(-c2  x  158. 5/T)  +  5.  (5) 

These  equations,  together  with  an  expression  for  the  line 
shape  f (u) ,  are  sufficient  for  calculating  the  limb  radiance 
profile  spectrally  resolved  over  the  147  jim  and  63  Jim  lines. 
The  Voigt  line  shape  was  used  in  this  work,  and  the  effects  of 
both  pressure  and  temperature  were  included  in  calculating  the 
collisional  (Lorentz)  width.  For  the  Lorentz  halfwidth  at  p  = 

1  atm  and  T  =  296K  we  used  the  value  0.329  cm  1 .  Pressure  has 

an  extremely  small  effect  in  the  inversions  to  be  described 
because  the  spectral  radiances  that  are  involved  correspond  to 
tangent  heights  greater  than  90  km  and,  more  importantly,  they 
have  been  degraded  by  an  instrument  spectral  scanning  function 
(see  Section  4)  that  produces  approximately  the  same  effect  as 
collisonal  broadening  corresponding  to  a  fixed,  high  pressure. 
The  inversion  utilizes  spectral  radiances  out  to  ten  Doppler 


halfwidths  from  the  line  center,  but  even  here  the  degraded 
spectral  radiance  has  a  very  weak  dependence  on  pressure. 
Pressure  broadening  would  have  a  significant  effect  for 
instrument  resolving  powers  very  much  higher  than  were 
considered  in  the  present  study;  see  Bullitt,  et  al.,  (1985). 
For  the  purpose  of  calculating  the  collisional  width  we  used 
the  pressure  profile  of  the  U.S.  Standard  Atmosphere  (1976) . 
Since  we  ignore  the  small  effect  of  differences  between  this 
reference  pressure  profile  and  the  "actual”  pressure  profile, 
the  calculated  limb  spectral  radiance  depends  solely  on  the 
vertical  profiles  of  the  O-atom  density  and  the  temperature. 
Hereinafter  these  profiles  will  be  referred  to  as  [0]  and  T, 
both  functions  of  altitude  H.  Both  the  monochromatic  limb 
radiance  profile  and  the  spectrally  resolved  limb  radiance 
profile  (the  monochromatic  radiance  degraded  to  the  finite 
resolving  power  of  the  instrument)  will  be  denoted  NV(HT), 
where  HT  is  the  tangent  height,  i.e.,  the  minimum  altitude 
along  the  line  of  sight.  We  will  identify  a  spectral  radiance 
or  other  spectral  quantity  as  monochromatic  or  spectrally 
resolved  where  necessary. 

For  later  reference  we  define  the  monochromatic  volume 
emission  rate, 

^([0], T)  -  A-[0]-f  (D-V0)gu  exp(-c2Eu/T)/Q(T)  ,  (6) 

which  is  the  number  of  photons  emitted  (in  all  directions)  per 
unit  time  per  unit  wavenumber  interval  by  a  unit  volume  whose 


O-atom  density  and  temperature  are  [0],T.  This  rate  equals 
471/ Az  times  the  last  term  of  equation  (1)  evaluated  in  the 
limit  as  Az  approaches  zero. 

The  altitude  variation  of  [0]  and  T  from  90  to  300  km 
altitude  is  typified  by  the  two  models  shown  in  Fig.  2. 

Figure  3  shows  the  monochromatic  limb  radiance  profile  NU(HT) 
computed  for  the  147  }im  line  using  atmospheric  model  ATOI4  of 
Fig.  2.  The  spectrum  for  HT  =  90  km  is  shown  as  a  dashed  curve 

for  greater  clarity.  Note  that  the  spectral  radiances  are 

,  2 
plotted  in  the  form  log  N^  vs.  (D-D0)  ,  which  results  in  a 

straight  line  when  has  the  Doppler  shape.  Over  the  range  of 

abscissa  in  Fig.  3,  which  corresponds  to  distances  \)-b0  less 

than  roughly  2.5  Doppler  halfwidths,  f (b)  is  well-approximated 

by  the  Doppler  shape.  For  the  largest  tangent  heights  Nu  has 

very  nearly  the  Doppler  shape  (same  shape  as  f (u) )  because  a) 

the  tangent  path  is  optically  thin,  and  b)  the  temperature  is 

nearly  uniform  over  the  path.  The  tangent  path  is  definitely 

not  optically  thin  below  HT  os  170  km. 

The  retrieval  technique,  basically  an  "onion-peeling” 
method  (Goldman  and  Saunders,  1979),  represents  the  atmosphere 
as  a  series  of  concentric  shells,  or  layers;  [0]  and  T  are 
treated  as  constant  within  each  layer.  The  retrieval  proceeds 
downward,  layer-by-layer,  from  some  maximum  altitude.  For  the 
purpose  of  discussion  the  sought-after  solution  may  be  regarded 
as  the  set  of  [0]s  and  Ts  for  each  of  the  layers  below  this 
maximum  altitude,  say,  300  km.  The  maximum  altitude  may  be  the 
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largest  tangent  height  HT  for  which  NW(HT)  has  been  measured, 
or  perhaps  the  largest  HT  for  which  the  signal-to-noise  of  the 
measurement  is  deemed  adequate.  Since  the  retrieval  technique 
involves  a  comparison  of  measured  and  calculated  limb  spectral 
radiances,  it  is  desirable  to  include  in  the  calculated  values 
the  contributions  of  layers  above  the  maximum  altitude.  This 
is  accomplished  by  a  starting  procedure  prior  to  the  onion¬ 
peeling  process. 

The  starting  procedure  exploits  the  fact  that  tangent 
paths  above  HT  ^  250  km  are  optically  thin,  and  that  NU(HT) 
tends  to  fall  off  exponentially  above  this  tangent  height  --  a 
consequence  of  the  nearly  constant  temperature  and  the 
exponential  decrease  in  [0]  with  altitude.  These  properties 
make  it  possible  to  analytically  invert  an  extrapolated  limb 
radiance  profile  via  the  Abel  transform  (see  Sharma,  et  al . , 
1988)  .  The  extrapolated  N^&j)  profile  corresponding  to 
Ht  >  300  km,  is  based  on  the  measured  NV(HT)  between  260  and 
300  km.  Specifically,  we  obtained  a  closed-form  approximation 
that  gives  the  spectrally  resolved  volume  emission  rate  at  H  = 
300  km  in  terms  of  Ny(HT  =  300  km)  and  a  single  radiance  scale 
height  derived  from  the  extrapolated  NV(HT)  profile.  This 
formula  is  used  to  recover  a  300-km  emission  rate  that 
corresponds  to  the  same  spectral  resolution  as  the  measured 
NU(HT)  .  '.en,  by  using  a  deconvolution  procedure  involving  the 

known  instrument  spectral  scanning  function,  the  monochromatic 
rate  ^ (K  -  300  km)  is  found.  A  corresponding  predicted  volume 
emission  rate  ^  (  [ 0 ]  ,  T )  can  be  obtained  from  equation  (6)  for 
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given  [0],T.  By  minimizing  the  least-squares  mismatch  between 
the  two  monochromatic  rates,  we  determine  [0]  and  T  at 
300  km.  These  two  values  define  [0}  and  T  at  all  higher 
altitudes  since  T  is  assumed  constant  above  300  km,  and 
because  [O]  vs.  H  is  assumed  exponential;  under  these 
conditions  and  for  the  optically  thin  case,  [0]  has  the  same 
scale  height  as  NU(HT)  . 

The  onion-peeling  procedure  is  pictured  in  Fig.  4  for  a 
representative  layer  at  some  tangent  height  less  than  the  300- 
km  maximum  altitude.  The  starting  procedure  and  previous 
onion-peeling  cycles  have  provided  [O]  and  T  for  all  layers 
above  the  one  in  question.  If  values  of  [O]  and  T  are 
assigned  to  the  layer,  one  can  calculate  the  limb  spectral 
radiance  for  the  layer  tangent  height.  The  solution  values  of 
[O]  and  T  are  defined  as  those  for  which  the  rms  difference 
between  the  calculated  spectrum  and  the  corresponding 
measured  spectrum  is  a  minimum.  No  assumption  is  made 
regarding  the  optical  thicknesses  of  the  individual  or 
aggregate  peeled  layers. 

The  physical  basis  of  the  retrieval  method  can  be  defined 
with  the  help  of  equation  (6)  for  ^([0],T).  The 
(monochromatic)  spectral  radiance  received  by  a  limb-viewing 
sensor  from  only  the  layer  at  the  tangent  height  is 
proportional  to  ^ ,  evaluated  for  the  layer's  [0]  and  T,  times 
the  spectral  transmittance  of  the  path  between  the  layer  and 
the  sensor.  The  received  single-layer  contribution  is  thus 
proportional  to  [O]  at  all  frequencies  v.  It  depends  on 
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temperature  through  the  line  shape  function  f  Cu)  and  the 
Boltzmann  factor  (normalized  by  Q(T)).  But  since  a  change  in 
the  latter  due  to  a  change  in  temperature  is  independent  of 
spectral  position,  the  retrieval  method  relies  on  f  (t)) ,  that 
is,  on  Doppler  broadening  of  the  line,  to  distinguish  between 
changes  in  temperature  and  O-atom  density.  Note  also  that  f  (U) 
decreases  with  temperature  near  the  line  center,  i.e.,  it 
opposes  the  change  in  the  Boltzmann  factor;  throughout  the 
wings  it  reinforces  the  Boltzmann  factor.  This  means  that 
there  is  a  range  of  temperatures  for  which  information  on 
temperature  is  obtained  mostly  from  the  line  wings.  Any 
Doppler  shift  due  to  platform  motion  and/or  thermospheric  wind 
along  the  line  of  sight  would  be  evident  in  the  spectrally 
calibrated  data  if  it  included  both  sides  of  the  line.  The 
shift  is  needed  to  correct  the  data  and  could  also  be  used  to 
determine  a  component  of  the  wind. 

Implicit  in  the  onion-peeling  nature  of  the  described 
inversion  method  is  the  assumption  of  a  spherically  stratified 
atmosphere.  The  tomographic,  two-dimensional  limb 
reconstruction  technique  described  by  Solomon,  et  al.,  (1984) 
could  be  applied  to  the  problem  were  it  not  for  the  fact  that 
reabsorption  is  significant;  we  are  not  aware  of  any 
multidimensional  limb  reconstruction  methods  without  the 
opt ically-thin  restriction.  A  one-dimensional  technique  that 
may,  in  the  present  application,  have  an  advantage  over  onion¬ 
peeling  with  respect  to  noise  propagation  is  briefly  described 
in  the  Conclusions  section. 


3.  ALGORITHMIC  DETAILS 


Let  Nm  (CALC)  and  Nvi  (MEAS )  denote  the  calculated  and 
measured  limb  radiance  spectra  at  wavenumber  ui  for  a 
particular  tangent  height  Ej.  Included  in  the  calculated 
spectrum  is  the  effect  of  an  instrument  line  shape  (ILS) /  or 
spectral  scanning  function,  representing  the  finite  spectral 
resolution  of  the  measurement.  Minimizing  the  sum-square 

y^[Nvl(CALC)  -Nui(MEAS)]2  (7) 


gives  [0]  and  T  at  altitude  H  =  HT .  We  elected  to  formally 
minimize  the  above  sum,  i.e.,  to  develop  coding  for  two 
function  procedures,  X { [ 0 ] , T }  and  Y{[0],T},  corresponding  to 
the  nonlinear  least-square  normal  equations,  symbolically, 

X  =  0  and  Y  =  0.  In  the  course  of  the  inversions,  these 
equations  are  solved  using  Newton's  method. 

The  atmosphere  must  be  finely  divided  to  obtain  reasonable 
accuracy  in  the  calculated  radiances  --  we  used  a  vertical 
thickness  of  one  km  for  the  layers.  Generally,  reducing  the 
layer  thickness,  which  reduces  each  layer's  contribution  to  the 
limb  radiance,  would  tend  to  increase  the  potential  for 
instability  in  the  retrieval  solution  due  to  noise/error 
propagation  if  the  solution  consists  of  independent  [0],T 
values  for  each  and  every  layer.  We  chose  to  solve  for  the  [0] 
and  T  of  just  certain  layers  spaced  more  than  one  km  apart, 
and  impose  the  constraint  that  log[0]  and  T  vary  linearly 
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between  these  one-km  "solution  layers".  The  distance  AH 
between  solution  layers  is  assigned  one  of  several  values 
depending  on  the  altitude  (see  Section  4  below).  The  [0]/T 
values  obtained  for  the  maximum  altitude  in  the  starting 
procedure  are  used  as  guesses  for  the  first  peeled  layer; 
thereafter,  guesses  are  obtained  by  linearly  extrapolating 
log[0]  and  T  from  the  previous  pair  of  solution  layers.  In 
the  event  the  normal  equations  have  no  solution,  the  guesses 
are  systematically  revised. 

The  profile  consisting  of  connected  temperature  solutions 
is  smoothed  above  K  =  130  km,  where  the  assigned  AH  is  10  km. 
Smoothing  is  performed  in  each  onion  peeling  cycle,  but  only 
for  altitudes  addressed  by  the  two  most  recent  cycles.  Thus, 
each  10-km  interval  above  H  =  130  km  is  smoothed  twice.  The 
smoothing  operation  eliminates  slope  discontinuities  and 
reduces  propagated  noise,  both  of  which  can  adversely  affect 
stability . 

4.  CHARACTERISTICS  OF  THE  SOLUTIONS 

Retrievals  were  performed  for  four  cases:  In  any  one 
case,  the  synthetic  spectral  radiance  data  represents  either 
the  63  )bn  or  147  nm  line,  and  atmospheric  model  ATOI2  or  ATOI4 
(Fig.  2) .  The  resolving  power,  R  =  Ug/Au,  where  Au  is  the 
full-width-at-half-maximum  of  the  ILS  used  to  degrade  che 
computed  monochromatic  radiances,  is  used  as  the  measure  of 
instrument  resolution.  We  used  an  ILS  corresponding  to  a  high 
finesse  Fabry-Perot  system.  Over  the  range  of  the  ILS  used  fo 


spectral  scanning,  it  is  well  approximated  by  the  Cauchy 
function  (Lorentz  shape) .  The  noise  equivalent  spectral 
radiance  (NESR)  ,  in  units  of  photons  s_1cm'2sr_1cm,  is  the  rms 
of  the  zero-mean  Gaussian  noise  added  to  the  spectrally 
degraded  radiances.  The  parameters  R  and  NESR  were 
independently  varied  over  a  wide  range  for  each  of  the  four 
cases.  The  following  discussion  emphasizes  results  obtained 
using  the  147  }im  line  since  those  obtained  from  the  stronger 
63  nm  line  were  generally  inferior  below  140  km  altitude. 

The  synthetic  147  Jim  data  for  each  tangent  height  consists 
of  25  samples  of  the  degraded  limb  spectral  radiance  over  half 
the  line.  The  spacing  of  the  samples  is  8  x  10”5  cm”1.  The 
half-width  of  the  infinitely  resolved  limb  spectral  radiance  at 
the  higher  tangent  heights,  for  model  ATOI4,  is  roughly 
2  x  10  4  cm  1  which  is  equivalent  to  2.5  samples.  Note  that  an 
instrument  half-width  equal  to  this  value  corresponds  to 
resolving  power  68/(4  x  10”4)  =  1.7  x  105.  Figure  5  shows  (for 
atmospheric  model  ATOI4)  the  monochromatic  147  |im  limb  spectral 
radiance  and  corresponding  degraded  spectral  radiance  for  two 
different  tangent  heights.  The  two  short-dashed  curves  are  Nv 
for  Ht=  90  km.  The  narrower  of  these  is  the  monochromatic  Nu, 
while  the  broader  one  is  the  degraded  Nw  for  resolving  power 
5  x  105.  The  solid  curves  are  similarly  defined,  except 
Ht  =  290  km.  The  long-dashed  curve  is  the  ILS  that  was 
convolved  with  the  monochromatic  Nws  to  obtain  the  degraded 
NyS .  The  points  in  Fig.  5  are  samples  of  the  degraded  Ny  for 
Ht=  290  km  after  addition  of  noise  corresponding  to  NESR  = 
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2  x  1011  ph  s_1cnf  2sr_1cm.  The  spectrally  degraded  data  for 
the  higher  tangent  heights  can  be  said  to  have  the  Voigt  shape, 
since  the  monochromatic  Nv  is  Voigt  (in  fact,  nearly  Doppler) 
and  the  ILS  is  approximately  Lorentz.  As  the  figure  indicates, 
the  degraded  Nu  typically  has  the  long  tail  of  the  ILS. 
Twenty-five  samples  accommodated  the  spectral  shapes 
corresponding  to  the  wide  range  of  resolving  powers  considered. 
In  practice,  a  sample  spacing  appropriate  to  the  actual 
instrument  resolving  power  would  be  used,  and  fewer  samples 
would  suffice. 

The  parameter  AH,  defined  earlier,  is  effectively  the 
vertical  resolution  of  the  sought-after  solution.  Initially, 
retrievals  were  performed  using  noise-free,  fully  resolved 
(R  =  oo)  data  and  AH  =  1  km  for  all  altitudes.  The  retrieval 
procedure  became  very  unstable  below  135  )cm.  At  approximately 
130  km  altitude  it  failed,  i.e.,  there  was  no  solution  to  the 
normal  equations.  This  behavior,  which  is  explained  below, 
gives  some  insight  into  the  effects  of  noise.  The  final 
version  of  the  retrieval  algorithm  uses  AH  =  10  km  above 
130  km,  AH  =  4  km  between  130  and  110  km,  and  AH  =  2  km  below 
110  km.  These  values  are  adequate  to  reproduce  the  true 
structure  of  the  model  [0]  and  T  profiles;  they  resulted  in 
solution  profiles  nearly  identical  to  the  model  profiles  over 
the  full  range  of  the  limb  radiance  data  (90  to  300  km)  for 
the  case  R  =  »,  NESR  =  0. 

The  spacing  in  tangent  height  of  the  data  used  in  the 
retrievals  is  the  same  as  the  altitude  spacing  AH  of  the 
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solution  layers.  The  results  reported  here  correspond  to  an 
instantaneous  field  of  view  (IFOV)  subtending  roughly  one  or 
two  km  at  the  tangent  altitude;  that  is,  IFOV  smearing  effects 
were  not  considered.  Synthetic  data  for  tangent  heights 
greater  than  220  km  was  generated  for  two-km  FLj.  intervals  and 
then  smoothed  to  obtain  the  spectral  radiances  at  10-km 
intervals.  A  given  NESR  refers  to  the  noise  level  before 
smoothing. 

The  solutions  of  the  nonlinear  normal  equations 
inevitably  contain  some  error,  even  if  the  data  has  no  noise. 
Errors  in  [0],T  propagate  downward  as  these  values  are  used  to 
compute  the  limb  spectral  radiance  Nu  for  a  lower  tangent 
height  and  then  to  retrieve  [0]  and  T  for  the  corresponding 
altitude.  Also,  as  the  retrieval  proceeds  downward,  the 
actual  temperature  at  the  tangent  altitude  changes,  and  the 
sensitivity  of  Ny  to  this  temperature  is  spectrally 
redistributed.  The  contribution  to  Ny  of  a  unit  volume  at  the 
tangent  altitude  depends  primarily  on  the  volume's  spectrally 
resolved  emission  rate  as  discussed  earlier.  Figure  6 
shows  versus  T  for  v  =  68  cm”1  (the  center  of  the  147  ^m 
line) ,  and  the  spectrally  nonresolved  rate  £  versus  T  for 
the  147  nm  line.  The  competing  temperature  dependencies  of  the 
Boltzmann  factor  and  the  lineshape  function  at  the  line  center 
result  in  the  broad  relative  maximum  in  £68(T)  near  T  =  500K. 

For  atmospheric  model  AT0I4  this  temperature  occurs  at 
approximately  130  km  altitude.  Thus,  for  HT=  130  km 
provides  temperature  information  about  altitude  H  =  130  km 


principally  in  the  wings  of  Nu.  But  a  significant  portion  of 
the  wing  radiance  is  due  to  higher  layers,  which  are  warmer 
(have  broader  due  to  Doppler  broadening)  than  the  130-km 
layer.  When  the  onion-peeling  procedure  has  reached 
Ej  —  130  km,  errors  in  the  computed  wing  radiance  have  become 
significant  in  comparision  to  the  actual  wing  radiance  from  the 
130-km  layer,  resulting  in  the  failure  mentioned  above. 
Obviously,  noise  in  the  data  can  be  similarly  propagated. 

It  was  found  that  the  solutions  tend  to  fall  into  three 
categories:  For  relatively  low  resolving  power  and  high  noise, 
the  retrieval  procedure  fails  at  approximately  130  km,  but 
provides  an  accurate  [0]  solution  and  a  usable,  albeit  noisey, 
temperature  solution  above  this  altitude.  Over  a  rather  wide 
range  of  higher  R  and  lower  NESR,  the  procedure  fails  at 
approximately  104  km,  but  yields  accurate  [0]  and  T  values 
for  H  >  104  km.  This  is  the  altitude  at  which  the  0-atom 
density  of  the  models  is  a  maximum.  Still  higher  R  and/or 
lower  NESR  results  in  very  accurate  retrieved  profiles  for  the 
entire  altitude  range,  90  to  300  km,  of  the  synthetic  data. 

The  failure  of  solutions  at  the  altitude  of  the  [ 0 ] — 
profile  peak  is  due  to  noise  propagation  in  combination  with 
high  atmospheric  opacity.  The  transmittance  from  the  tangent 
point  to  the  exoatmospheric  observing  instrument,  at  the 
center  of  the  147  ^m  line,  is  less  than  ^0.01  for  tangent 
heights  less  than  104  km.  The  63  |im  line  gave  poorer  results 
than  the  147  line,  particularly  for  the  lower  altitudes, 
because  of  its  h.iaher  opacity.  The  63-|im  line-center 
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transmittance  to  the  tangent  point  is  less  than  ~0.01  for 
tangent  heights  less  than  120  km. 

5.  SUMMARY  OF  RESULTS  (for  147  Line) 

The  noise  levels  and  resolving  powers  corresponding  to  the 
three  solution  categories  are  defined  in  Fig.  7.  For  example, 
any  combination  of  R  and  NESR  corresponding  to  a  point 
between  the  two  lowest  solid  curves  results  in  accurate  [0] 
and  T  solutions  between  300  km  altitude  and  the  peak  of  the 
[0]  profile,  near  104  km.  Figure  7  is  a  composite  of  two  such 
plots  obtained  separately  for  the  two  modal  atmospheres.  Thus, 
it  may  not  accurately  characterize  the  solution  when  the  true 
[0]  or  T  profile  is  very  different  from  those  of  the  models. 

The  dashed  curves  in  Fig.  7  are  lines  of  constant  maximum 
signal-to-noise,  defined  as  the  maximum  of  the  measured  Ny(HT) 
for  u  =  uQ  =  68  cm’1,  divided  by  the  NESR.  This  is  also  the 
required  dynamic  range  of  the  measurement  system.  The  maximum 
signal  occurs  at  HT  120  km  for  a  given  resolving  power.  As 
R  is  increased  the  maximum  signal  becomes  larger,  and  a  given 
maximum  S/N  corresponds  to  larger  NESR.  The  signal-to-noise  is 
considerably  lower  than  the  maximum  in  the  wings  of  the  line 
and  at  the  higher  tangent  heights,  i.e.,  Nu (HT)  has  a  very 
large  range  of  values.  Figure  5  shows  that  for  NESR  = 

2  x  1011  ph  s’1cm"2sr’1cm  and  R  =  5  x  105,  the  lowest  S/N  in 
the  data  is  unity,  and  that  the  NU(HT)  data  spans  almost  four 
decades . 

Figures  8  through  10  show  retrieved  [0],T  profiles 
representative  of  the  three  categories.  The  point  symbols  in 


Fig.  7  identify  the  NESRs  for  Figs.  8-10/  the  resolving  power 
is  5  x  105.  Figures  8  and  10  correspond  to  model  ATOI2, 
representing  a  cool  thermosphere;  Fig.  9  corresponds  to  model 
ATOI4 . 

Figures  9  and  10  show  clearly  that  the  [0]  profile  can  be 
recovered  to  high  relative  accuracy  even  when  there  are  large 
errors,  e.g.,  100K,  in  the  retrieved  temperature  profile. 

This  is  true  at  least  for  the  higher  altitudes.  Generally, 
the  temperature  errors  decrease  from  high  to  low  altitude,  due 
to  the  increase  in  S/N,  provided  the  solution  doesn't  fail. 

When  solutions  can  be  obtained  down  to  90  km  (see  Fig.  8),  the 
temperature  errors  at  the  lower  altitudes  (90  to  ~104  km) 
become  very  small,  but  the  errors  in  [0]  tend  to  be  larger 
below  the  (0]-atom  density  peak  than  above  the  peak. 

6 .  CONCLUSIONS 

The  described  technique  will  yield  the  translational 
temperature  and  0- atom  density  profiles  from  a  limb  radiance 
scan  spectrally  resolved  over  the  147  line.  The  level  of 
noise  in  the  data  and  the  spectral  resolution  have,  of  course, 
an  effect  on  the  accuracy  of  the  retrieved  profiles,  but  more 
importantly,  they  determine  the  altitude  at  which  the  onion-peel 
retrieval  technique  becomes  unstable.  Very  accurate  profiles 
can  be  recovered  from  ~  300  km  down  to  90  km  altitude  if,  for 
example,  the  resolving  power  is  5  x  105  and  the  noise- 
equivalent  spectral  radiance  is  less  than  7  x  1010 
ph  s"1cm"2sr“1cm.  With  half  this  resolving  power  and  ten  times 


the  noise  it  is  possible  to  recover  an  accurate  O-atom  density 
profile  and  a  rough  approximation  to  the  temperature  profile 
above  ~  130  k,n  altitude.  Above  130  km  the  instrument  need  only 
resolve  approximately  10  km  at  the  tangent  altitude,  whereas 
two-km  vertical  resolution  is  required  to  retrieve  the  sharp 
peak  of  the  O-atom  density  profile  near  100  km  altitude.  The 
63  ^m  line,  which  is  more  opaque  than  the  147  line,  yields 
inferior  solutions  below  ~  140  km  altitude. 

A  resolving  power  of  5  x  105  and  NESR  of  7  x  1010 
ph  s'1cm'2sr~1cm  at  wavelength  147  pjn  work  out  to  a  required 
noise-equivalent  radiance  of  2  x  10  14  W  cm~2sr  1 .  These 
signal  requirements,  and  the  corresponding  two-km  vertical 
resolution,  were  translated  to  basic  design  requirements  for  a 
Fabry-Perot  spectrometer  system  operating  at  orbital  altitude. 
It  was  found  that  the  etaion  mesh  reflectors,  assuming  a 
parallel  mesh  design,  would  have  to  be  much  larger  in  size  than 
we  believe  is  practicable.  However,  a  confocal  design  based  on 
spherical  mesh  reflectors,  which  was  suggested  for  the  present 
application  by  H.  A.  Smith  (personal  communication,  1987),  can 
theoretically  provide  the  required  performance  with  an  etaion 
measuring  a  few  centimeters.  While  the  etaion  size  is  more 
reasonable  in  the  confocal  design,  it  remains  to  be 
demonstrated  that  the  theoretical  resolving  power  of  5  x  105 
can  be  achieved  in  practice.  With  the  confocal  design,  and  a 
stressed  Ge:Ga  detector  providing  a  noise  equivalent  power  of 
5  x  10'U  W  Hz’1/2,  and  a  one-Hz  noise  bandwidth,  the  system 
would  need  cold  foreoptics  approximately  a  meter  in  diameter. 


The  one-Hz  noise  bandwidth  corresponds  to  the  measurement  of 
Nu ( H? )  at  a  single  combination  of  u  and  HT.  The  time 
available  to  scan  u  and  HT  may  make  it  necessary  to  use 
multiple  detectors  and/or  a  larger  collector.  Thus,  retrievals 
down  to  90  km  by  the  described  technique  imply  a  somewhat 
advanced  sensor.  A  similar  system  with  larger  IF OV  and  cold 
foreoptics  of  ~15  cm  diameter  would  allow  retrievals  from  300 
km  down  to  130  km  altitude. 

It  is  likely  that  the  reported  signal-to-noise  and/or 
spectral  resolution  requirements  could  be  reduced  somewhat  by 
replacing  the  onion-peeling  procedure  by  a  method  that  fits, 
simultaneously,  the  measured  limb  spectral  radiances  for  many 
tangent  heights  in  the  data  set.  This  is  aptly  called  the 
global-fit  approach  by  Carlotti  (1988) .  Chang  and  Weinreb 
(1985)  refer  to  it  as  the  Levenberg-Marquardt  method  after  the 
algorithm  they  used  to  effect  the  nonlinear  least-squares  fit. 
Both  papers  compare  results  obtained  by  the  global-fit  and 
onion-peeling  methods,  and  conclude  that  the  global  technique, 
which  does  not  tend  to  propagate  noise,  yields  better 
solutions.  The  paper  by  Carlotti  describes  a  retrieval  of 
stratospheric  ozone  concentrations  by  the  global-fit  approach 
that  required  substantial  execution  time  on  a  Cray  computer. 
This  discouraged  us  from  applying  the  global-fit  method  to  the 
present  problem,  which  we  estimate  would  require  even  greater 
computer  resources.  Also,  we  wanted  to  first  compare  the 
spectrally  resolved  and  non-resolved  approaches  for  temperature 
and  O-atom  density  retrieval,  based  on  onion-peeling.  As 


pointed  out  by  Chang  and  Weinreb  (1985),  the  onion-peeling 
method  is  valuable  as  a  diagnostic  rool;  in  the  present 
application  it  helped  us  recognize  a  fundamental  sensitivity 
problem  near  130  km  altitude  (T  —  500K) .  Onion-peeling 
combined  with  global  fitting,  where  the  latter  is  used  only 
for  a  selected  range  of  altitudes,  e.g.,  90  to  130  km,  could 
conceivably  be  used  in  our  application  to  achieve  reduced 
sensitivity  to  noise  together  with  reasonable  computing  time. 
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FIGURE  CAPTIONS 

1. )  Approximate  energy-level  diagram  for  the  atomic  oxygen  01 

transitions.  Values  are  given  for  the  energies  and  for  the 
transition  lifetimes  and  Einstein  A-coef ficents  (after 
Fischer  and  Saha,  1983) . 

2. )  Two  models  for  the  vertical  profiles  of  atomic  oxygen 

density  and  translational  temperature.  They  represent  a 
cool  thermosphere  (AT0I2)  and  a  relatively  warm 
thermosphere  (AT0I4) . 

3. )  The  spectrally  resolved  monochromatic  limb  radiance  computed 

for  various  tangent  heights  for  the  147  nm  line  and 
atmospheric  model  AT0I4.  The  abscissa  is  proportional  to 
the  square  of  distance  from  the  line  center. 

4. )  The  geometry  of  the  onion-peel  retrieval  technique. 

5. )  Monochromatic  and  spectrally  degraded  limb  radiance  for 

tangent  height  90  km  (short-dashed  curves)  and  tangent  height 
290  km  (solid  curves) .  The  degraded  spectral  radiances 
correspond  to  resolving  power  R  =  5  x  105.  The 
instrument  line  shape  for  this  R  (long-dashed  curve)  is 
normalized  to  fit  the  radiance  scale.  The  spectral 
radiances  represent  the  147  line  and  atmospheric  model 
ATOI4 . 
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6. )  The  volume  emission  rate  evaluated  at  the  center  of 

the  147  jam  line  (u  =  68  cm  l)  ,  and  4(T)  for  the  147  jim 
line.  Both  are  normalized  to  their  values  at  T  =  296K. 
The  spectrally  resolved  rate  ^  was  calculated  from 
equation  (6),  in  which  the  Doppler  shape  was  used  for 
f  (V-\>0)  . 

7. )  Noise-equivalent  photon  spectral  radiances  (NESR)  and 

resolving  powers  at  X  =  147  Jim  corresponding  to  the  three 
solution  categories.  The  dashed  curves  indicate  the 
corresponding  maximum  signal-to-noise  (required  dynamic 
range)  . 

8. )  Actual  and  retrieved  [0]  and  T  profiles  for  resolving 

power  5  x  10".  The  NESR  is  7  x  1010  ph  s_1cm”2sr_1cm. 

3.)  Same  as  Fig.  8,  except  the  NESR  is  2  x  1011 

-1  -2  -1 
ph  s  cm  sr  cm. 

10.)  Same  as  Figs.  8  and  9,  except  the  NESR  is  5  x  10U 

-1  -2  -1 
ph  s  cm  sr  cm. 
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